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The Eastern Diamondback Rattlesnake, Crotalus adaman¬ 
teus Palisot de Beauvois, 1799, is the largest rattlesnake 
species native to North America, reaching an average adult 
length of 130 cm (Conant & Collins 1991). The species 
is distributed throughout the southeastern United States, 
although populations are mostly located in Florida, Ala¬ 
bama, and southern Georgia. Despite a significant popu¬ 
lation decline since the 1950s (Timmerman & Martin 
2003), mostly as a result of habitat loss, habitat fragmenta¬ 
tion, and indiscriminate killings by humans (e.g., “rattle¬ 
snake roundups”, Means 2009), the species has not been 
afforded specific protection. Additionally, limited data are 
available on the most fundamental aspects of C. adaman¬ 
teus biology, including population biology, life history, and 
population genetics (Hoss et al. 2010). As urbanization 
and other anthropogenic effects continue to degrade and 
destroy viable habitat, small and isolated biological pre¬ 
serves may be the last resort for this and similar species. 
However, snakes with large home ranges (such as C. ada¬ 
manteus) may not benefit from small preserves because 
they are unlikely to meet the requirements for stable popu¬ 
lations (Dodd et al. 1988, Timmerman 1995) and natural 
population genetic mechanisms may be compromised. For 
example, the population on our university campus, Florida 
Gulf Coast University (FGCU, Fort Myers, Florida, USA), 
may be rather small and isolated because of surrounding 
development. If so, inbreeding could lead to lowered het¬ 
erozygosity and fitness. 

Despite the need for conservation of Eastern Diamond- 
backed Rattlesnakes, few genetic studies have been per¬ 
formed. For example, only two polymorphic microsatellite 
markers have been identified for C. adamanteus (Bushar 


et al. 2001) and none have been developed specifically for 
this species. Because of successful cross-species amplifi¬ 
cation of PCR primers in Crotalus and other snake taxa 
(Bushar et al. 2001, Stapley et al. 2005, Anderson 2006, 
Lukoschek & Avise 2011, Pozarowski et al. 2012), we 
aimed to test published microsatellite loci from four con¬ 
generic rattlesnakes with C. adamanteus DNA and charac¬ 
terize the successful loci in the FGCU population. Estab¬ 
lishing polymorphic molecular markers in this species will 
allow advanced studies in behavioral ecology and conser¬ 
vation genetics. 

We extracted adult genomic DNA from scale clips, 
blood, or liver samples with standard phenol-chloroform 
procedures. We initially tested one DNA sample for cross¬ 
amplification success with 25 published microsatellite 
primer pairs from four other rattlesnake species (C. atrox 
and C. scutulatus, Pozarowski et al. 2012; C. tigris, Gold¬ 
berg et al. 2003, Munguia-Vega et al. 2009; and C. viridis, 
Oyler-McCance et al. 2005). Oligonucleotide primers 
were labelled for fluorescent detection using a three-prim¬ 
er reaction system as in Boutin-Ganache et al. (2001) 
and Croshaw & Glenn (2003). The forward primer of 
each pair had a 5’ M13 tag added (s’-GGAAACAGCTAT- 
GACCATG-3’), and a third dye-labeled primer (6FAM) 
was included in each reaction. We used standard touch¬ 
down polymerase chain reaction (PCR) (Don et al. 1991) 
with the following reaction conditions: 25.0 pL total reac¬ 
tion volume with 12.5 uL GoTaq Master Mix (2x, Promega), 
0.16 pM forward primer, 0.16 pM reverse primer, 0.2 pM 
6FAM-labeled M13 primer, and 20-50 ng DNA template. 
Based on visualization of PCR products with 2% agarose 
gels, we selected 10 loci for further study. Reactions were 
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Table 1. Characterization of seven microsatellite markers from N = 15 individual Crotalus adamanteus in southwest Florida, USA. 
Each forward primer had the M13 tag (5’-GGAAACAGCTATGACCATG-3’) appended to the 5’ end. The sample size (n), number of 
alleles observed (K), expected and observed product sizes (bp), TouchDown temperature (TD) in thermal regime, observed heterozy¬ 
gosity (H o ), expected heterozygosity (H e ), and Nei’s unbiased gene diversity (H s ) are also listed for each locus. P is the p-value based 
on the results of a statistical test for Hardy-Weinberg equilibrium in the FGCU population (Genepop). Expected product size comes 
from the original source publication. Locus Ca274 was monomorphic and therefore offered no information for statistical analyses. 


Locus Source Species Repeat motif Primer sequence (5’-3’) n K Expected Observed P TD H H H 

product product 
size (bp) size (bp) 


CS234 Pozarowski C. scutulatus 
et al. 2012 


Cal22 Pozarowski C. atrox 
et al. 2012 


Ca281 Pozarowski C. atrox 
et al. 2012 


Ca238 Pozarowski C. atrox 
et al. 2012 


Crti23 Munguia- C. tigris 
Vega et al. 


(GAAA) 5 (GGAA ) 4 F: CAGATGGTTAAAAG- 
(GAAA) 20 GGAA ggcatca 


(GAAA ) 4 


R: GGAATTCAACCAT- 
AAAAGT GACA 


(TG ) 15 F: CCCCCAGGAAA- 

GAATATGTATG 


14 10 438-546 455-508 1.000 55 0.929 0.886 0.919 


15 2 187-225 189-197 0.983 60 0.667 0.515 0.533 


R: CACATCTTCAGGTT- 
GCCAAA 

(GAAA) 19 TAAA F: CCCAGCACACT- 12 7 253-455 322-355 0.372 60 0.583 0.786 0.820 

(GAAA ) 4 CACAATGTC 

R: GCTCCTAGCCAGGAAT- 
GTCA 


(AGGC) n (AGGA ) 7 F: CTCTGCTCTGCCATT- 15 6 315-411 356-384 0.599 60 0.867 0.818 0.846 
TCACA 

R: GAGCCAGTGCTTT- 
GTTTTCC 

(GT) 12 (GA ) 12 F: GATTGTGTGGTGTT- 15 2 224-234 199-203 0.060 60 0.133 0.515 0.533 
TATTGTTGC 


2009 


Cal31 Pozarowski C. atrox (TC ) 19 
et al. 2012 


Ca274 Pozarowski C. atrox (GT ) 21 
et al. 2012 


R: GGATGCTCCATCACTA- 
GACG 

F: GGGATTCAG- 14 9 187-340 227-275 0.426 55 0.714 0.767 0.795 

GCCAATTTTTA 

R: TGGTAATTTTGCATT- 

TCAGCA 

F: AGCCGCTTGTCT- 15 1 292-334 295 N/A 65 0 0 0 

CAAATTGT 

R: TGTCCTCTGATTT- 

GCCCATT 


then run on polyacrylamide gels (Mini-Protean Precast, 
Biorad) to visually screen for polymorphism. Because 
three primer pairs did not show consistent amplification in 
the expected size range, we ran the samples from just seven 
loci on an ABI 3130 automated DNA sequencer (Univer¬ 
sity of Florida’s Interdisciplinary Center for Biotechnology 
Research) with LIZ 600 as an internal size standard. We 
called alleles with Peak Scanner vi.o (Applied Biosystems). 

After confirming the presence of scorable peak patterns, 
we then optimized thermal and reaction conditions before 
genotyping 15 total snakes collected in southwest Florida, 
USA [11 adults and one offspring from the Florida Gulf 
Coast University campus (FGCU) two adults from Florida 
Panther National Wildlife Refuge, and one adult from Ko- 
reshan State Park]. Reactions were run with three different 
touchdown PCR thermal cycling programs on T100 Ther¬ 
mal Cyclers (Biorad). The amplification conditions consist¬ 
ed of programs using a range of annealing temperatures 
from either 65-55°C, 6o-50°C or 55-45°C. Each program 
consists of an initial denaturation of 96°C for 3 min, fol¬ 
lowed by 5 cycles of 96°C for 20 s, the highest annealing 


temperature for 30 s, and extension of 72°C for 1 min, fol¬ 
lowed by 21 cycles of 96°C for 30 s, the highest annealing 
temperature minus o.5°C each cycle for 30 s, and 72°C for 
1 min, and finally 10 cycles of 96°C for 30 s, the lowest an¬ 
nealing temperature for 30 s, and 72°C for 1 min. 

Finally, the complete genotypic data set (N = 15 snakes) 
was used to characterize the loci in Cervus 2.0 (Marshall 
et al. 1998 ). We tested each locus for Hardy-Weinberg equi¬ 
librium (HWE, exact test) and linkage disequilibrium in 
the FGCU population (N = 11 adult snakes) with Genepop 
V 4 . 1.4 (Raymond & Rousset 1995 ) using all default set¬ 
tings. 

Of the seven consistently amplifying microsatellite loci, 
six were polymorphic (Table 1). For polymorphic loci, 2 to 
10 alleles (mean = 6.0) were amplified. Observed hetero¬ 
zygosity ranged from 0.13 to 0.93 (mean = 0.65) and Nei’s 
unbiased gene diversity ranged from 0.53 to 0.92 (mean = 
0.74). HWE was not rejected in any case (p-value ranged 
from 0.06 to 1.00) but one locus in particular (Crti23) 
showed high heterozygote deficiency and had a high es¬ 
timated frequency of null alleles (0.58). The null hypoth- 
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esis of independence of genotypes could not be rejected for 
any locus pair (p-value ranged from 0.22 to 1.00). The en¬ 
tire marker panel showed low non-exclusion probabilities 
(0.015 for the second parent and 0.000003 for identity). 

Our results are consistent with the relatively high cross¬ 
species amplification success reported within the Crotalus 
genus and in other snakes (King 2009, Lukoschek & 
Avise 2011, Rice et al. 2017). Most of the 25 primer pairs we 
tested generated clear products that appeared to be in the 
expected size range on agarose gels. Further assessment of 
the applicability of existing primers should be performed 
for additional species within the Crotalinae and other 
groups to provide markers useful in population genetics. 
Moreover, with six of the seven markers we assayed being 
polymorphic in C. adamanteus, it is likely that establish¬ 
ing panels of robust, informative microsatellite loci will not 
be especially difficult in this species. Additional advanced 
research applications in population genetics, conservation 
genetics, and behavioral ecology of this species, which are 
sorely lacking, are possible with these markers. However, 
one locus (Crti23) had substantial heterozygote deficiency, 
likely because of null alleles. This primer pair should be 
tested in other populations and with larger samples to de¬ 
termine its utility. Nevertheless, our study offers the possi¬ 
bility for neutral evolutionary processes to be added to the 
extensive work on the adaptive evolution of venom pro¬ 
teins in this species (e.g., Margres et al. 2014, 2016). 

Despite our concern that the small campus population 
of C. adamanteus at FGCU could have low genetic varia¬ 
tion, our results revealed levels of gene diversity (mean = 
0.74 for FGCU snakes) and polymorphism (mean = 4.5 al¬ 
leles per locus) that were consistent with other population 
genetic studies of Crotalus (King 2009). Perhaps this pop¬ 
ulation is either larger than we know or is connected to 
healthy populations on the east side of the FGCU campus, 
where development is less extensive. Alternatively, if pop¬ 
ulation declines were recent they might not yet have re¬ 
sulted in detectable losses of genetic variation. Much more 
data are needed before robust population genetic conclu¬ 
sions can be made, and further ongoing analysis may re¬ 
veal answers to these questions. The microsatellite DNA 
loci reported here allow for future studies concerning con¬ 
servation and population genetics, paternity analysis, and 
individual relatedness within wild C. adamanteus popula¬ 
tions. 
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